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f"^ ■ Abstract 

1 We present new results for the matrix elements of the Qe and Q4 penguin operators, evaluated 

f"S |' in a large- N c approach which incorporates important 0{N%-J-) unfactorized contributions. Our 

approach shows analytic matching between short- and long-distance scale dependences within 
dimensional renormalization schemes, such as MS. Numerically, we find that there is a large 
pH I positive contribution to the A J = 1/2 matrix element of Qe and hence to the direct CP- violation 

parameter e'/e. We also present results for the A7 = 1/2 rule in K — * nn amplitudes, which 
incorporate the related and important "eye-diagram" contributions of 0(N^^-) from the Q2 
, operator (i.e. the penguin-like contraction). The results lead to an enhancement of the AI = 1/2 

« I ■ effective coupling. The origin of the large unfactorized contributions which we find is discussed 

' in terms of the relevant scales of the problem. 



Oh: 



1 Introduction 



The so called penguin operators in Particle Physics have a long history. Their existence was first 
pointed out by the ITEP group yQ who showed that in the process of integrating out heavy fields 
in the electroweak theory of the Standard Model, in the presence of QCD interactions, there appear 
new local four-quark operators, like e.g. (summation over quark colour indices within brackets is 
understood): 

<?4 = 4 (*L'l ll qL)(qLlvd L ) and Q 6 = -8 Y (sLto)fe^), (1-1) 

besides the two conventional operators 

Q2 = 4(sl7 m ul)(ul7 m c?l) and Q 1 = 4(s I y i d L )(u L 7 M w L ) , (1.2) 



which had been considered up to then [21 EJ - These four-quark operators, modulated by their ap- 
propriate Wilson coefficients Cj (/i) , which are functions of the masses of the fields which have been 
integrated out and the scale [i of whatever renormalization scheme has been used to carry out this 
integration, are part of the effective Hamiltonian which describes the weak interactions of quarks at 
intermediate energies of a few GeV 1 , 

„ 10 

H cS (AS = 1) = -J=V* a Vud dMQi(n) • (1-3) 

The ITEP group also claimed that matrix elements of the Qq penguin operator were particularly 
important and could be at the origin of the experimentally well established AT =1/2 enhancement in 
K — > 7T7T transitions. However, the dynamical mechanism they proposed was shown to violate Current 
Algebra Ward identities El CI IE] • In fact, in the approximation where the Qq operator factorizes 
into a product of bilinear quark densities, the K — > h-k matrix elements turn out to be too small to 
explain the AT =1/2 enhancement. 

In the meantime, it was also shown that the Wilson coefficient Ce(/i) acquires a large imaginary 
part, as a result of the integration of the heavy t and b quarks and the flavour mixing structure of 
the CKM matrix in the Standard Model. As a result, the evaluation of K — > irir matrix elements of 
the Qq operator, besides its potential contribution to AI =1/2 transitions, has also become of utmost 
importance, as an ingredient for a plausible understanding of the observed size of the CP-violation 
parameter e'/e in the Standard Model. 

The effective Lagrangian which describes \A.S\ = 1 transitions, like K — > 7T7t and K — > 7T7r7r, in 
the presence of electromagnetic interactions and to lowest order in chiral perturbation theory (%PT) , 
O(p ) and 0(p 2 ) in this case, has the following structure : 

£ cff =1 = ~% VM W C& + g ^ £27 + g2g ^ tr (u^UQr)] , (1.4) 

where 

£ 8 = Y, ( £ m)2* (^)ia and £ 27 = \{C^) 21 (£") 13 + (£^23 (^)n • (1.5) 



with 
and 



i=l,2,3 



4 32) = S i3 S j2 , Q L = Qr = Q = diag.(2/3, -1/3, -1/3). (1.7) 

The pion decay coupling constant Fq is the one in the chiral limit where the quark masses u, d, s 
are neglected (Fq ~ 87 MeV). The matrix field U collects the Goldstone fields of the spontaneously 

1 See e.g. the lectures of A. Buras in ref. Q], where earlier references can also be found. 



£ p = -zF 2 U{x)^D^U{x) , (1.6) 
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broken chiral symmetry of the QCD Lagrangian with three massless flavours, and D^U denotes the 
covariant derivative: D^U = d^U — ir^U + iUl^, in the presence of external chiral sources and r M of 
left- and right-handed currents. Notice that the octet term proportional to g 8 induces pure AI =1/2 
transitions, while the term proportional to g 2 7 induces both AI = 1/2 and AI = 3/2 transitions. The 
coupling constants g 8 , g 27 and g ew encode the dynamics of the integrated degrees of freedom in the 
chiral limit. These include, not only the heavy quark flavours, but also the hadronic light flavour 
states (resonances) other than the Goldstone fields explicitly present in the xPT Lagrangian 

The coupling constant g 8 has both a real part (relevant for the AI =1/2 rule) and an imaginary 
part (which induces direct CP violation) . The main purpose of this work is to present an evaluation of 
the contributions to both e'/e and AI =1/2 induced by the QCD penguin operators in Eq. (|l.l(l . as 
well as from the so called eye-like configurations of the Qi operator in Eq. I|1.2f) . We shall do this within 
the framework of the 1/N C expansion. The methodology 2 is the same as the one which has been applied 
to other calculations of similar low-energy observables reported elsewhere |^ El El EI El E] • It 
allows for important improvements with respect to earlier calculations within the framework of the 
1/N C expansion, reported in ref. Q2J (at leading order) and in refs. |18l 1191 l2TH 1211 122| (at next-to- 
leading order). 

The paper is organized as follows. Section 2 collects the general formula from which the bosoniza- 
tions of the penguin operators can be obtained. Section 3 gives the results for the usual factorized 
contributions. Section 4 contains a detailed discussion of the calculation of the unfactorized contri- 
butions of Q{N c rif). Section 5 contains the corresponding analytic results for the coupling constant 
g g with both the factorized and unfactorized contributions. The phenomenological implications for 
the AI =1/2 rule and for e'/e are presented in sections 6 and 7 respectively. Finally in section 8 we 
present a discussion of our results and conclusions. 



2 Bosonization of the Penguin Operators 

The Dirac operator T> x of the QCD Lagrangian 

C QCD = -±G$G^ + iqV x q, 
in the presence of external chiral sources r M , M. and is defined as follows: 



1 - 75 



1 + 75 



1 1 + 75 

2 



(2.1) 



(2.2) 



where G M denotes the gluon matrix field, and G$ the gluon field strength tensor (a = 1, . . . , 2V| — 1). 
The bosonization of the penguin operators Q4 and Qq, to 0(N%) and 0(N C ), is then formally defined 
by the functional integrals [231 



(Q 4 (x)) =4Trp- 1 H 
+4 



x 

74, 



x — TV V~ x {—%\ — 



V J (27r)4 e { l) Sl,(x) 3j - 



6V y m . SV X 



(2.3) 



and 



(Qe(x)) = -8 Tr - - Tr V~\-i 



x 



* v '5Ml(x) 3j x v '5M{y) o2 



+8 



2 See e.g. refs. |9lllf)| for recent 



4y f d q c -iQ-(x-y) T ( X)-i(-i) 5 -2a x>- x f i) 5T>X \ (2 4) 
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where the trace Tr here also includes the functional integration over the planar gluonic configurations 
which leave the quarks at the edge 0] The first terms in Eqs. (|2.3J) and (|2.4(l correspond 

to the factorized contributions, illustrated in Fig. 1, which are O(N^); while the second terms in 
Eqs. H2.3J) and l|2.4|) correspond to the unfactorized 0(N C ) contributions, also illustrated in Fig. f . 
These unfactorized contributions involve integrals over the incoming q momenta, which are regulariza- 
tion dependent. For consistency, they have to be defined in the same renormalization scheme as those 
of the corresponding Wilson coefficients. With p the conjugate momentum operator, in the absence 
of the external chiral sources, the full quark propagator in x space is given by the expression 



where 



1> + T [Ic^+r^] 
d 



M 



1-75 



JWti±ig 



\y) 



(x\ j\y) = 7^ 



dxt 1 



G 



S(x - y) . 



(2.5) 



(2.6) 



The bosonization of four-quark operators in xPT is then obtained via an appropriate chiral expansion 
in powers of the l a , r a , M. and external sources in the propagators. 






factorized 



unfactorized 



Fig. 1 Factorized and unfactorized contributions corresponding to Eqs. \2.!-i\) and j2.4\ ). The crosses 
correspond to the explicit variation of the Dirac propagator with respect to external sources ( terms like 
&M(y) 2 * n d3P i> ^ e s °tid ti ne s represent quarks, in the presence of soft external chiral sources, 
propagating in a planar gluon background. 



Let us recall that the spectrum of hadronic states which can contribute to these Green's functions 
in QCD in the large-7V c limit |241 12"5] consists of an infinite number of narrow stable meson states 
which are flavour nonets. In the real world, however, the physical hadronic states have widths, 
which are subleading in 1/N C . The fact that the ratios of widths to masses are of 0(1/N C ) is at the 
very basis of the phenomenological success of the large-iV c approach. Indirectly, this also suggests 
that in calculating physical observables, the contributions which stem from narrow states are the 
ones which may be potentially important. Notice that in the case of the bosonization of four-quark 
operators, the large- N c QCD spectrum of narrow states contributes to the factorized as well as to the 
unfactorized patterns. Consequently one might already expect the unfactorized piece to be sizeable. 
Our calculation shows that this expectation is confirmed and the unfactorized contributions from these 
penguin operators turn out to be larger than those from the factorized terms. We shall later argue 
that this result becomes rather natural when formulated in terms of the different hadronic scales 
involved. 



3 The Factorized Contributions 

Let us first discuss the factorized contributions from the Qg operator. The operators (slIr) and 
(to^L) are Noether densities of the QCD Lagrangian in Eqs. (|2.1|l and (|2.2|l : 
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I- \ -I \\ 1+ ^ 5 ( \ -n ( \ ^QCdO) , ql s 
(SLqR) = q(x)X 3j — 2~Q\ X ) = D sg{x) = - 5M f^ 3 . i 3 - 1 ) 

/ — 7 \ -( 1-75 , X ^ / \ ^QCD( X ) ,„ 

(<7L£«fl) = q{x)X j 2 — — q{x) = Gq d {x) = — 5M ^ x y 2 ■ ( 3 - 2 ) 

On the other hand, the QCD effective Lagrangian which describes the strong interactions of Goldstone 
particles at low energies is given, in SU(3)l x SU(3)r, by a string of terms 

C x = ^F 2 {trtDpUDfrf) + 2Btr(MU^ + UM ] )} 
+2BL 5 tr[D^U' i D> 1 U(M' i U + U ] M)] + • ■ • 

+L 3 tr(Z? M C/ 1 D"U D„U j D U U) + iL 9 tr(Fg v D^UD^ + F^ V D V U^D V U) + ■■■ , (3.3) 

where the first line gives the well known 0(p 2 ) terms, and only those terms of the 0(p 4 ) Lagrangian |77] 
which will be needed in this paper have been explicitly written out. Notice that we are considering 
a framework in which the massive rj' particle has been integrated out; while, strictly speaking, in 
a large-iV c formulation one should consider a U(3)l x U(3)r description of the chiral Lagrangian 
with an explicit ?7o _ singlet Goldstone particle. We have checked that for the observables discussed 
in this paper, and at the level of the 0(rif/N c ) approximation that we are retaining, there is no 
phcnomenological difference between the two formulations. 

The bosonization of the Noether densities D Sq (x) and D q d(x) can then be readily obtained as 
variations of the chiral Lagrangian in Eq. (|3.3|l with respect to the same external sources: 

Ds " {x) = ~ SM*(xl ■ ^ WtTXsj { \ F % U W + WD^DOU +•••}+ 0(p 4 ) , (3.4) 
G M = - //r*/ 3 ? 2BtrX j2 {- A F*U\x) + L 5 D^D»Utf + ■■■}+ 0(p A ) , (3.5) 



6M (x) j2 



where, among the 0{p 2 ) contributions, we have only written the term induced by the L5 coupling of the 
0(p 4 ) chiral Lagrangian because it is the only term which contributes to the factorized bosonization 
of the Qq operator. Indeed, although the two densities D Sq (x) and D q d{x) start with terms of O(p ), 
their product -because of the fact that the U matrix is unitary- only starts at 0(p 2 ) in the chiral 
expansion and, because of the flavour structure, this product can only depend on the L5 coupling. 

On the other hand, the factorized contribution from the Q4 operator is rather straightforward. 
With 

/- n \ -, a 1 - 75 / \ 7-p / \ ^QCP(^) ,„ 

(sLj^qL) = q[x)X3jT—^—q(x) = L? q (x) = - , (3.6) 

(- \ -( w u!-75 M _ rfl m ^QCd(») ,„ 
(qLTd L ) = q{x)X j2 T—^—q{x) = L^ d {x) = - — , (3.7) 

the bosonization of these left-current densities, to lowest order in the chiral expansion starts at 0(p) 
and are given by the expressions: 

L%{x) => \F 2 tr [X 3j (D^)U] + 0(p 3 ) (3.8) 
LU X ) MD»tf)U] + 0( P 3 ) . (3.9) 
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The contribution to the coupling constant g 8 from the factorized patterns of the Q4 and Qq 
penguin operators can now be readily obtained from the definitions in Eqs. (|1.3|1 . (|1.4|l and the results 
in Eqs. I|3.4j) . (|3.5|) and (|3.8|) . I|H.9|) . One thus obtains the well known result 3 



I factorizcd 

i,Q4 



c 6 (m) 



-16L, 



W) 2 1 



F 2 F 4 



C 4 (m)1 • 



(3.10) 



where we have used the fact that 



B=i2l(#)l 



(3.11) 



It is well known that in this approximation, polychromatic penguins don't fly, as emphasized in ref. [8], 



4 Beyond Factorization 

Because of the chiral structure of the lowest order electroweak Lagrangian in Eq. I|1.4|) , we have to 
expand the Dirac propagators in the unfactorized patterns in Eqs. (j2.3(l and Ij2.4[l up to two chiral 
powers in the external sources. There is, however, a certain arbitrariness in the choice of the external 
sources to expand. The simplest (and recommended) choice is the one which makes it explicit that all 
the perturbative QCD (pQCD) short-distance contributions are factored out in the Wilson coefficients 
Ci(/z). Since the operators Q4 and Qq transform like 8 L x 1 R operators under chiral rotations, 
expanding the propagators in right-handed sources will necessarily bring in Green's functions with 
extra i? M = g7 M 1 ~ l 2 75 q currents only, and there is no way then that the product of these (1 L x 8 fi ) 
operators with the initial (8 L x 1 R ) operator can produce a 1 L x 1 R term, which ensures no 
mixing with possible pQCD contributions in the long-distance evaluation. Technically, this can be 
accomplished by an appropriate "tout a droite" rotation 4 of the quarks fields so that in the rotated 
basis the right-handed quark field has absorbed all the Goldstone degrees of freedom and the Dirac 
operator has as a chiral vector connection the external field only. One way to implement this is as 
follows: define Ql and Qr quark fields, such that 

Ql(x) = £, L (x)q L {x) and Qr{x) = £ R (x)q R (x) , (4.1) 

where £l(x) and £r(x) are left and right SU(3)l x SU(3)r/ SU(3)v coset representatives , with 
£r( x )£l( x ) — U(x). Then, choose a gauge where 

fo(a!) = l, while t R (x)=U(x). (4.2) 

The Dirac operator in this "tout a droite" rotated basis is rather simple 

D x = r (d M + zGV - ilJ—^\ + U^D^Ur^ 1 + W X M + iM*U . (4.3) 

The chiral expansion, in the chiral limit where M. — A4' = 0, can then be made in powers of the 
external source which collects the full Goldstone structure; i.e., the term U^D^U"/^ "t 75 in Eq. (|4.3J) . 

4-A Determination of the unfactorized Green functions 

As already mentioned, we are only considering in this work the bosonization of the unfactorized terms 
in Eqs. H2.3JI and Ij2.4|l of 0(N 2 ^) in the 1/N C expansion, and in the chiral limit. As we shall sec, 
the Physics features which already emerge at that level deserve attention. It is possible, however, to 
extend the calculation in the chiral limit to the rest of the 0(N 2 jj-) contributions, which are not 
enhanced by a nj factor, something which we plan to do in the near future. 

3 See e.g. refs. |8lll7ll28l where earlier references can also be found. 

4 In the french cycling jargon, "tout a droite" means to put the front and rear gears so as to reach maximum speed 
(the biggest front plateau with the smallest rear socket). This is what we recommend to do here with the Dirac operator! 
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Fig. 2 "Eye"-like configurations contributing to the unfactorized bosonization of the Qq operator. 
The solid dots in the quark propagator correspond to external right-handed sources of soft quark 
currents. 



Let us discuss the case of the Q§ operator in some detail. The only configurations which can con- 
tribute in the approximation that we are considering, are the so called "eye" -like configurations [2"F)| . 
where the two right-handed sources come out from the same propagator in Eq. (|2.4(l , as illustrated in 
Fig. 2. They produce the functional integral 5 

d4 1 f j4 . m-/_l 1 + 7s i 



(Q 6 )U = 8 n f J J d*y Tr(x\ 2 '° j lv)* 

which gives a term with the required chiral structure: SnfXtvX 3 2U^D a UU^DfjU = —infXtvX^DcJJ^DpU, 
modulated by the integral (recall Eq. (|2.6[l ) 

ll. [ d\Tv(x\ l+l5 1 \v)(v\ ±—£ 1 7° 1 + 75 1 ^ 1+75 1 \x) (4 5) 

The integral in question can be related to a well defined QCD Green's function, which is the connected 
Green's function 

W DG f R (q,k) = }im i 3 f d 4 xd 4 yd 4 ze^e^ye- lk - z (0\T{D Sq (x)G q MRt(y)RL^ (4-6) 
where (no summation over the flavour index q in Eq. (|4.6|l ) 

D s<1 (x)=q(x)X 3j ^^q(x), G 9<J (0) = g(0)A J2 ^*g(0) , (4.7) 

and 

RUV) = ^V)^il al -^q{y) , RL{z) = q(z)X 13 ^^^q(z) . (4.8) 

This Green's function can be viewed as a two-point function of the two D Sq (x) and Gqd(0) operators 
which define Qe, with two soft insertions of the composite operators R^ u (y) and R^ s (z). The relation 
to the integral in Eq. 1)4. 5J1 . and hence to (Q6)|c y c, is the following 



(Q 6 )U. = -tr\ 32 D Q U^D U x 8n/ x U — -^-1 



(4.9) 



5 The contribution from the term where the first propagator in Eq. 12.41 is expanded instead of the second one, is 
modulated by the flavour factor ^ ^j2^Zj = 0, and therefore it vanishes. 



G 



with the integral restricted to the term proportional to g a " . In fact, this term appears naturally, when 
performing the integral over the solid angle dfl q , which has the form (Q 2 = —q 2 ): 

J dn q W DG f R (q, k) = - W DGRR {Q 2 ) , (4.10) 

where the transversality in the four-vector k follows from Current Algebra Ward identities. We are still 
left with an integral of the invariant function Wdgrr(Q 2 ) over the full cuclidean range (0 < Q 2 < oo) 
which has to be done in the same renormalization scheme which has been adopted when doing the 
calculation of the short-distance Wilson coefficient Cq(p) in Eq. Ijl.^J) . i-e. in the MS scheme. 

The case of the unfactorized bosonization of the Q4 operator is rather similar. To 0(p 2 ) it is 
governed by a connected Green's function which can be viewed as a two-point function of the two 
L^ q {x) and L^ d (0) operators which define Q4, with two soft insertions of the composite operators 

R du(y) and R as(z); i-e., 

(Q 4 )U = ttX^Dal^DpU x Anj x [g^ J ^W^(q,k) 

where here 

W L ^(9,fc) = |im z 3 | d 4 *Ad 4 ^-V^e-^(0|T{L^^ (4.12) 

with no summation over the flavour q index in the currents. The Ward identities again ensure that 
the integral over the solid angle in Eq. 14.11f) also depends on one invariant function only, 

g, u J dQ q W£Z£g{q, k) = - g°A W LLRR (Q 2 ) , (4.13) 

and the term proportional to g a/3 appears then naturally, as well. 

It is now possible to express the contribution to the coupling constant g% of the electroweak chiral 
Lagrangian in Eq. I|1.4|l from the factorized patterns of the Q4 and Qq penguin operators to 0(N 2 ), 
as well as from their unfactorized patterns of 0(N 2 ^-), which appear as integrals of the invariant 
functions Wdgrr{Q 2 ) and ~Wll RR {Q 2 ), in the following way: 



(4.11) 



jo,.,,,. = C 6 (/i) <j-16L 5 — — + 



(4V)£/2 r dQ 2 (Q 2 y-^ 2 W DGRR (Q 2 ) 



+ C 4 (/i) 1 



An 



f 



16ir 2 F$ 



{A-KfJL 



2\e/2 



r(2 - e/2) J 



r(2-e/2)y 

dQ 2 {Q 2 ) 1 ^' 2 W LLRR {Q 2 ) 



MS- 



MS 



(4.14) 



Notice that the integrals over Q 2 are MS renormalized; i.e., poles in 1/e as well as the constant 
log47r — are removed. Let us remember, however, that the scale dependence of the renormalized 
constant L${n) in xPT is 



L s Qj,) = L 5 (M P ) - 



167T- 



•log 



JL 



with r,5 = 



(4.15) 



and coupling constants in xPT have been renormalized in the Bern MS scheme where poles in 1/e, 
as well as the constant log47r — 75 + 1, are removed. This difference of MS renormalization schemes 
has to be taken into account consistently in order to exhibit the overall p scale cancellation between 
the short- and long-distance contributions, at the order of approximation that we are working in the 
1/N C expansion. Here we have chosen to renormalize all the UV divergences with a single scale /z e . 
Therefore, in Eq. I|4.14|l there is a /Lt-scale dependence both in parameters such as L 5 ([i) and (-ipifj) (/j,) , 
as well as in the short-distance Wilson coefficients, such as Cq^(/j,). There is, however, freedom to 
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rcnormalize L§ with a different scale, which could be called ^ chiral as conventionally done in xPT. 
That this is possible is due to the fact that the renormalization of L5 does not mix with that of 
the quark-condensate, or the Wilson coefficient. Notice also that, in the MS renormalization scheme 
which we are using, where quark masses (and hence the pseudoscalar masses) are set to zero, there 
are no chiral loop contributions to the factorized term in Eq. (|4. 14|l . The fact that all the UV scale fi 
dependence cancels out in the end in the coupling constant g% is a nontrivial check that the interplay 
between short- and long-distances has been brought under control. 



4-B Why one could expect large unfactorized contributions 

At this level, we think it necessary to have a short discussion explaining the relative importance of 
the factorized versus unfactorized contributions in Eq. I|4.14|) . It has long been recognized that the 
calculation of the coupling constant g g is a multi-scale problem. Already at short distances, the 
perturbative running of the Wilson coefficients has to deal with several different scales from M\y 
down to the charm mass. However, it is perhaps not so much emphasized that, even in the long- 
distance regime, there are also two distinctive scales playing a role in Eq. (|4. 14|l . Although in QCD 
all scales are ultimately related to A^r, it is a fact of life that Fo is much smaller than a typical 
hadronic scale, of the size of an average resonance mass, Mr ~ 1 GeV. Also the typical scale defined 
by a QCD vacuum condensate 6 , A c , is known to be smaller than Mr. One then has that, roughly, 
F -C A c <C Mr. Because Fq is of 0(N C ) while M R is of 0(1) a naive interpretation of the large-A^ 
counting could lead one to conclude that Fq ^> M R , which is obviously wrong. It is this hierarchy 
of scales which is at the origin of why the unfactorized contribution, although of Q(N C ), may be as 
large as -or sometimes larger than- the factorized one, even though the latter is formally of 0(N 2 ). 
For instance, we already found that in the case of the calculation of Bk in the chiral limit |12l I16| . 
the unfactorized contributions are of the order of 50%. In the present case of the QqjQh penguin 
operators, the unfactorized contributions actually turn out to be much larger than the factorized ones, 
as we shall soon see. We conclude that in the case of four-quark operators, a naive interpretation of 
the 1 /N c expansion may be misleading and has to be done after the contributions from the different 
physical scales have been separated. 

In the case we are considering, knowing that (as will be confirmed in the next section by explicit 
calculation): 

/>oo />oo 

/ dQ 2 Q 2 W DG rr(Q 2 ) ~ B 2 F 2 , and / dQ 2 Q 2 W L lrr(Q 2 ) ~ F 2 M 2 , (4.16) 
Jo Jo 

up to logarithmic factors, and that -L5 ~ Fq/M r , the ratio of the unfactorized to the factorized 
contributions in Eq. 14.14J) is given by precisely the ratio of scales 



unfactorized M 



factorized 16tt 2 Fq 



(4.17) 



which is a number of 0(1), although large- N c suppressed. This is true for both the contributions from 
Q 6 and Q 4 to g$ in Eq. (|4.14l) . 

Is this a shortcoming of the large— N c expansion? Not necessarily. It is rather that the standard 
analysis, as used, e.g., in ref. |25j was only applied to situations in which the two scales Fq and Mr 
do not compete for the control of the large- A*^ expansion, as it is the case in Eq. I|4.14|l . In fact, we 
think that there is nothing specially new about this dynamical effect of scales. Even the perturbative 
expansion in QED can be affected by the presence of two widely separated masses, because they give 
rise to large logarithms of their ratio. In the case of the electroweak theory, Bardeen et al. in ref. [2j 
were the first to point out that the large- N c expansion is also affected by the presence of the large value 
of log-Mv^/m-charm- The point that we would like to emphasize here is that, also in the low-energy 
region, ratios of scales such as M R / (4ttFq) 2 may upset a large- N c counting done in a naive way. On 
the contrary, once the contributions from these two scales have been clearly separated, it should be 
safer applying a large- Af c argument. For instance, one might fear that since the unfactorized term is 

6 We are thinking here of something like (■tpip) 1 ^ at fi = 1 GeV. 
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larger than the factorized one, the whole expansion breaks down and the next subleading terms will 
be even bigger. Although a definitive answer to this point can only come from the solution to QCD in 
the large--/V c expansion, which is not known, we would like to argue that there is no reason to expect 
this breakdown. Subleading terms to those of Eq. Ij4.14|l are expected to give contributions of the 
order of the width of the resonances involved, or of the violation of the Zweig rule. There is presently 
no evidence against these effects being reasonably small except, perhaps, in some very specific cases 
where the effect of the strong Goldstone-Goldstone interaction in the J = I = channel appears, 
which may deserve special attention and which we plan to study elsewhere. 



4-C Long and short distance constraints on the functions Wdgrr(Q 2 ) 

AND W LLRR (Q 2 ) 

The low-Q 2 behaviour of these functions is governed by x^T. An explicit calculation gives the results: 



m*W DGRR m x + (-L 5 + - i3 ) W + 0(Ct,), (4.18) 

\hnW LLRR (Q 2 ) = -0 + (-y£ 3 + ^ 9 ) + 0(Q 2 ). (4.19) 



The double pole contribution to Wdgrr{Q 2 ), the L$ contribution to Wdgrr{Q 2 ) and the simple pole 
contribution to V^llrr{Q 2 ) agree with the result first obtained in refs. |18M20| and ^5] respectively 7 . 
To our knowledge, the other terms have not been calculated before. 

On the other hand, the high-Q 2 behaviour of the functions YVdgrr{Q 2 ) and Wllrr(Q 2 ) is 
governed by the operator product expansion of the D &)G density currents in Eq. (|4.6J) and the 
L®L currents in Eq. (|4.12|) . with the results: 



ja.iwtf) - 4 ^-ir'?^ ,4 - 21) 



We remark that in Eqs. (|4.18|) - (|4.21|) the chiral couplings L5.9 are the ones corresponding to leading 
order at large- N c and, therefore, do not run with the scale. 

In this work we shall content ourselves with only the leading-log approximation. This allows us 
to do away with the evaluation of the terms of ©[e^r] in the OPE, which means that we shall not be 
able to specify scheme dependences in our MS-renormalization calculation of the long-distance effects. 

In large-iV c QCD, the Green's functions Wdgrr(Q 2 ) and Wllrr{Q 2 ) are meromorphic functions 
in the Q 2 variable. The most general structure they can have is: 



Q 2 Wdgrr(Q 2 



1 (BF ) 2 



4 \ "* a i , 6 \ "* A I 8 \ "* li (a 99 n 

Mhad. 2^Q2 + M 2 + Mhad. 2v <Q2 + M f)2 + ^had. 2v fQ2 + M 2)3 ' 

Q 2 W LLRR {Q 2 ) = Mhad. Yl Q2 + M 2 + ^had. (Q2 + M 2 ) 2 + ^ ad ' ^ (Q 2 + M 2 ) 3 ' 



where fihad. is an arbitrary hadronic mass scale, so as to make the residues on, Pi, ^i and oip /3'-, 
j'j dimensionless. This scale can be interpreted as the scale at which the four-quark Lagrangian 

7 The Lz contribution to Eq. 14.181 was omitted in ref. |20l for reasons which turn out to be incorrect. 
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in Eq. is matched onto the chiral Lagrangian of Eq. (|1.4fl . i-e. the scale at which the meson 

resonances are integrated out. Physical results, however, do not depend on the particular choice of 
the scale phud., except for higher order terms. The sums are, in principle, extended to an infinite 
number of states, as illustrated in Fig. 3. In practice, it is useful to start with the minimal hadronic 
approximation 8 (MHA) where the number of non-Goldstone states is limited to a minimum number 
of lowest mass states, with appropriate quantum numbers; the minimum number which is necessary 
to satisfy the leading OPE constraints. In our case, the MHA we shall adopt is the one with a 1~ 
vector pole of mass My and a + scalar pole of mass Ms- 




Fig. 3 Different types of hadronic tree diagrams which can contribute to the Green's functions 
Q 2 Wdgr.r(Q 2 ) Q 2 Wllrm{Q 2 ) (replacing x — * I) in the large-N c limit. The indices i, j, k label 
the possible poles. 



With the notation (M v = 770 MeV and M s ^ 1 GeV) 



Pv = ^~ and = (4.24) 



^had. ^had. A'had 

the MHA parameterizations are, 

8 VMhad./ z z + Pv z + p s 

and 

-<L«1(^Ld.) = + , (4-26) 

z + pv z + ps 

with the residues ay, as, ot' v and a' s solutions of the system of equations defined by the short-distance 
constraints 

-(BFq) +Mhad. 2^ Ul = ~q n ~ F o +3 77 — p2 (4.27) 



i=V,S 



j=V,S 6 6 



'See e.g. ref. 1101 and references therein for details 
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and the long-distance constraints 




1 2 3 4 1 2 3 4 



Fig. 4 The functions zYVdgrr [z[i 2 lad ) (with the Goldstone pole removed) and zVV ^™ RR \{z /^ ad ) in 
Eqs. and \4-.2b\j , versus z = Q 2 / [i\ ad for the particular choice fih&d. — 1 GeV and | (-0-0) l 1 ^ 3 — 

250 MeV. The values of the other input parameters are specified in the Appendix. Notice also the 
different vertical scales in the two plots. 



The shapes of the resulting functions in Eqs. (|4.25|l and (|4.26(1 are shown in Fig. 4. A few comments 
are now in order: 

• The reason why we are showing the shape of the functions zW^q RR {z) and zW L ^™ R K R {z) , instead 
of wj^ RR (z) and W^™ RR (z) themselves, is that these are precisely the relevant integrands which 
appear in the integrals in Eq. (|4.14|) . The pion pole in z x W^ RR (z) does not contribute to 
the "area" in dimensional regularization, and this is why we can remove it from the integrand 
with impunity. 

• In spite of the fact that the Goldstone pole (i.e. the term proportional to (-BFb) 2 in Eq. (|4.25|l ) 
does not contribute to the Q 2 integral in dimensional regularization, the sum of the residues of 
the hadronic simple poles on is constrained by the residue of this Goldstone pole through 
the OPE, as shown in Eq. (|4.27l) : and in fact, it is this term which largely dominates the short- 
distance constraint. This is the subtle way in which the lowest order chiral Lagrangian, which 
does not contribute at the factorized level, shows its presence at the unfactorized level. It is 
precisely this term which is responsible for the renormalization of L5. 

• The residue of the single Goldstone pole contribution to Wdgrr is fixed by a combination of 
couplings of the 0(p 4 ) chiral Lagrangian, as seen in Eq. (|4.18(l . The same residue provides a 
constraint on the sum J^i ~ °f the hadronic parameters, as seen in Eq. (|4.29|l . 
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• Numerically, the chiral factor L5 — \L 3 ~ 10 -2 which appears in Eqs. H4.18|l and (|4.29l) is rather 
large, and it is at the origin of the fact that the unfactorized contributions will turn out to be 
so important, since it is this quantity (times (tpip) 2 ) which fixes the value at the origin of the 
integrand Q 2 x W DGRR {Q 2 ) in Eq. jOl^ . 

• The corresponding chiral slope factor — ^§L 3 + |ig ~ 30 x 10~ 3 for Wllrr(Q 2 ) in Eq. (|4.19jl is 
even larger. However, in this case, the intercept at the origin is fixed by the 0(p 2 ) term — §Fo 
(see Eq. H4.30| l1. which is of a reasonable size. 

The MHA parameterizations of the functions WdGrr(Q 2 ) and yV^™ RR (Q 2 ) may, however, be 
improved. On the short-distance side, they only approach their OPE behaviour at rather large 
Q 2 values. On the long-distance side the behaviour of y^LLRfiiQ 2 ) at small Q 2 fails to reproduce 
the known slope in Eq. (|4.19(1 . This is why we also went beyond the MHA and considered more 
sophisticated parameterizations of the general structure given in Eqs. (|4.22|1 and H4.23[l . demanding 
that they interpolate smoothly between the known chiral behaviour and the OPE behaviour. In the 
case of Wdgrr {Q 2 ) we included a vector simple pole, a scalar simple and double poles, and an excited 
pseudoscalar pole. For Wllrr(Q 2 ) we included a scalar simple pole together with simple, double and 
triple pole for a vector, an excited vector, and an axial-vector states. The extra information to fix 
the residues of the poles was obtained by demanding that Wdgrr{Q 2 ) and WbLRRiQ 2 ) reproduce 
the OPE at various points, Q 2 > 9 GeV 2 . This scale was arbitrarily chosen as being low enough to 
help producing a smooth interpolation but high enough for the OPE behaviour to be trusted. Typical 
examples of the shapes we get for nfzWoGRRiQ 2 ) an d nfzWLLRR{Q 2 ) are shown, respectively, 
in Figs. 5 and 6 below (the thick solid curves). For the purpose of comparison, we also show in 
the same plots the OPE behaviour (the thick dashed curves), as well as the xPT slope in the case 
of zWllrr(Q 2 ) (the thin dashed line in Fig. 6). The slope of zWdgrr(Q 2 ) (with the pion pole 
removed) could be determined from the couplings of the 0(p 6 ) chiral Lagrangian. The calculation, 
which is rather involved, is beyond the scope of this paper. 




Fig. 5 Shape of zWDGRR(z/J> 2 la d) (pion pole removed) versus z (solid curve) for a large-N c type 
parameterization as in Eq. \4-Z 6 ^ which matches smoothly the short-distance OPE behaviour (the 
thick dashed line) with the long-distance xPT behaviour. The curves correspond to /ihad. = 1 GeV, 
^5|iar g c-]v c = 10~ 3 and I (?/;?/>) I 1 / 3 = 250 MeV; the values of the other input parameters are specified in 
the Appendix. 
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Fig. 6 Shape of zWllrrIz^^ ) versus z (solid curve) for a large-N c type parameterization as in 
Eq. which matches smoothly the short-distance OPE behaviour (the thick dashed line) with the 

long-distance X^T behaviour i.e., the point at the origin and the slope (the thin dashed line). The 
curves correspond to /ihad. = 1 GcV, L 5 |ia rB e-iv c = 1CT 3 and |(^)| 1/3 = 250 MeV; the values of the 
other input parameters are specified in the Appendix. 



5 Analytical results for the coupling constant g 8 from QCD Penguins 

We are now in a position to do the integrals in Eq. (|4.14|) and hence to obtain an evaluation of the 
contribution to g 8 from QCD penguins, beyond the factorization result in Eq. <|3.10|) . We shall keep 
the two-loop evaluation of the Wilson coefficients C4 and C§, as obtained e.g. in ref. down to 
the charm mass scale m c = 1.3 GeV; while the evolution from m c to an arbitrary hadronic scale 
/ihad. < TTt c is done at the leading log approximation, which is consistent with the bosonization of 
the Qi and Qq operators we have done. This results in the following expression for g 8 | q 6 ,q 4 valid to 
leading and next-to-leading order in the 1/N C expansion, including terms of O(jf-): 

s §\q 6 ,Q4 = C e(m c )f 6 (m c ; /ihad.) + C 4 (m c )/ 4 (m c ; /i ha d.) , (5.1) 

where 



fe(m c ; /ihad.) 



a s {m c ) 

a S (Mhad.) 



SI / a s (m c ) 



-8n f 



Mhad. 

16tt 2 F 4 



( a* log pi 



had J 

pi 



-16L 5 (// h ad. 



Ill 

2P? 



In/ 

9N r . 



XY y/ ^had. 

F 6 
1 

a s (m c ) 
a s (^had.) 



(5.2) 



and 



fi(m c ; /ihad.) 



q s (m e ) 

,a S (/ihad.) y 



I "/ 

II N c 



1 + An 



Phad. \ " 

167T 2 K 4 ^ 
u 3 



1 , % Hi 
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In/ 

9N n 



a s (m e ) 

a S (Mhad.) 



■16L. 



(5.3) 



It is worthwhile to compare the functions fe(mc] A*had.) and fi(fn c ; /ihad.) with the corresponding 
expressions obtained from factorization; i.e., 



/ 6 (m c ; /i ha d.) 



facto 



a s {m c ) 
a s (Mhad.! 



-16L 5 (/w.) ^ 



^had. 




/4(m c ;/i h ad.)L. = 1 



I facto 



(5.4) 
(5.5) 



We observe that f% (m c ; /ihad. ) | f act0 rizod ^ as a ra th er strong dependence on the choice of the hadronic 
scale /ihad. (a fact which is often ignored in phenomenological discussions on e'/e in the literature), 
while /e(Tn c ; //had.) has a remarkably smooth dependence. This behaviour is illustrated in Fig. 7, 
where we plot these functions, normalized to their respective values at /ihad. = 1 GeV; i.e., 

7 / \ ffti^Tlc ;/ihad.) i _? / \ /s (^cj /ihad. )factorized / r ~\ 

/ 6 (m c ;^had.) = 7-7 : - n and / 6 (m c ; /ihad.Jfactorized = 7-7 1 - 77: (5.6) 

/ 6 (m c ; 1 GeV) h{m c ; 1 GoV)f act orizod 

versus /ihad. in the range 0.8 GeV < /ihad. < 1.3 GeV. The reason why we normalize the /e functions 
to their value at /ihad. = 1 GeV is that, in absolute value, fe{m c ; /ihad.) turns out to be larger than 
/ 6 (w c ;/ihad.)factorizod; for example, at /ihad. = 0.8 GeV, we find that 

fe(m c ; 0.8 GeV) ~3x f 6 (m c ; 0.8 GeV) 

factorizcd j 

(5.7) 

which is an important enhancement 9 . We also find a similar enhancement, though perhaps less 
dramatic, of the unfactorized contribution to f^(m c ; /ihad.), which as we shall see later, is a welcome 
feature towards a phenomenological understanding of the observed AI = 1/2 rule. 



f 6 ( m c; ^had.) 



0.8 0.9 1 1.1 1.2 1.3 

^had. 

Fig. 7 The dependence on the choice of the /ihad. scale of the functions f6(m c ; /ihad.) (the flat solid 
curve) and f6{m c ; /ihadjfactorized (the dotted curve), normalized to their respective values at /ihad. = 
1 GeV. See the definitions in Eq. j5.6j) . 



9 This is in qualitative agreement with the numerical results found by the authors of refs. |T1 22 within the framework 
of the extended Nambu-Jona-Lasinio model (ENJL). 



1.75 
1.5 
1.25 
1 

0.75 
0.5 
0.25 



f fi ( m c; fhai), , . . 
«. « v ' factorized 
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6 Phenomenology of K — > nn Amplitudes and the AI = 1/2 Rule 

We shall define the decomposition of physical K — > irir amplitudes into isospin amplitudes Ai—op, as 
follows: 

A[K° -» n+TT-} = iA e tS ° + -^-=iA 2 e U2 , (6.1) 

A[K° -> 7r°7r°] = zAoe l<5 ° - V2iA2e Wa , (6.2) 
tt+tt ] = |iA a e Wa , (6.3) 

where, 5o and 52 are the J = irir phase shifts with 1 = and I = 2 at the if mass. The chiral 
Lagrangian in Eq. I|1.4|> gives the following contributions to the Aj isospin amplitudes 



(6.4) 



Gp T , T „ [5 ,, ,1 Ts 11 



^2 = 2F Q |^g27(M| - m 2) _ __ e 2 g ^| . (6 . 5) 

In the large-iV c limit, when only the factorized contributions are taken into account, there is a dy- 
namical symmetry which relates the g 8 and g 27 couplings: 

g 8.lo(AT2) = S27|e,(jv2) = g i ( 6 -6) 

while we recall that, experimentally, 

| gi U. ~ 5.1 and |g27|c XP . ~ 0.29 : (6.7) 
although the values to be explained can be reduced to 

|gsU. -3.3 and \g^\^. ~ 0.23 , (6.8) 

if one takes into account the enhancement already provided by the calculated 0(p 4 ) chiral correc- 
tions no no cm EE]. 

The unfactorized 0(N C ) contributions break the dynamical symmetry which is at the origin of 
the disastrous prediction in Eq. I|6.6|) : but, as was pointed out in ref. |32], there still remains a 
smaller dynamical symmetry at that level of approximation. This symmetry relates the weak matrix 
elements of the Q2 and Q\ four-quark operators in Eq. 1|1.2[) (neglecting their mixing with the penguin 
operators) to those of the AS* = 2 operator 

Qas=2 = (s L rd L )((s L ^d L ) , (^°|QAs= 2 (0)|iY O ) = f K M\ 3 as= 2 (m) , (6.9) 

in the following way: 

Rc S8|q 2 ,Q X = + ^5as=2(m)^) + Mti - j?SAS=2(/u)J , (6.10) 

and 

3 

S27 = + z 2^)] g5AS=2(^) , (6-11) 

where we are using the same definition of the Wilson coefficients as in Buras et al. i.e., 

C i (p)=z i (ri+Ty i {n), with T = -^±. (6.12) 

V us V u d 

Our present work updates this dynamical symmetry with three new ingredients: 
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1. The calculation of the constant <?as=2(m)j within the same framework of the l/N c expansion 
discussed here, which has been reported in refs. [121 116) . 

2. Our new result in Eq. (|5.1(l adds an extra contribution to the real part of the g 8 coupling 
constant 

Reg s |q 6 ,q 4 = z &{v)h{v ; Mhad.) + z A {p)f A {p, ; fi had .) , (6.13) 

3. The contribution to Reg 8 from the eye-like configuration of the Qi operator, at 0(N C ), can be 
read off straightforwardly from our calculation of the hadronization of the Q4 penguin operator, 
and adds an important extra contribution 

Reg 2.| Q2 , Eye = Z2 ^~ > ^ had -) ~ ^ • ( 6 ' 14 ) 

This eye contribution is essential since it provides the log// dependence which cancels with 
the ones from ze(fJ,)fe(fi; /ihad.) and Z4(//)/4(/x;/ihad.) 10 - Note that this calculation is to our 
knowledge the first one where such a scale dependence cancellation in the Qi~Q^fi mixing 
sector is explicitly shown. 

Altogether, and at the 0(N C ) we are working, we get 

R-egg = ZlO) ^-1 + ^5AS=2(/i)^ + Z2(H) ^1 - ^9AS=2(^ + Z ^)^j [AOSMhad.) - l] 

+^6(j")/6(M!A t h a d.) + Ziip) fi(fl 5 /"had.) ■ (6.15) 

The relation in Eq. I|6.11fl which fixes g 27 in terms of <7as=2(m) has been known for a long time |34| 
and, in the chiral limit, it holds to all orders in the l/N c expansion. In our case, the numerical result 
for the invariant Bk factor ^1 ^5] obtained in the chiral limit 

B K = jC A s= 2 (m) x .9as= 2 (m) = 0.36 ± 0.15 (6.16) 

implies 

g27 = 0.29 ±0.12, (6.17) 

which is perfectly consistent with the experimental value given in Eq. (|6.7|) . This is reassuring, 
because, at the order of approximations we are working, it is a complete calculation. However, for 
phenomenological applications, it remains to be seen how this result will be modified in the presence 
of chiral corrections. 

By contrast, the result for Reg g in Eq. H6.15JI . is not yet a full 0(N C ) calculation. We think, 
however, that it is worthwhile to present the numerical results which one already obtains at this level. 
There is no contribution to Reg 8 from the Q4.6 operators for [i > m c . Therefore, at fi = m c ~ 1.3 GeV, 
only the first line in Eq. (|6.15(l is nonvanishing. The corresponding numerical results we obtain for 
Reg 8 1 p , for two input values of | (^tp) (/x = 2 GeV) | (see discussion below) and letting the large-N c 

value of L5 vary in the range 1 x 10~ 3 < L^\i &tst ^ Nc < 2 x 10 -3 , have been tabulated in Table 1 below. 
For each entry, the range of the results is the one corresponding to Fig. 6 by varying /ihad. in the 
interval 0.8 GeV < /ihad. ^ 1.3 GeV, while, at the same time, allowing for violations of factorization 
in the residue of the OPE in Eq. (|4.21|) by an extra factor of 2. This extra factor should take into 
account the fact that we are finding large deviations from the factorization of matrix elements of 
four-quark operators and, therefore, there could be large corrections also in the residues of the OPE 
in Eq. (H23J. 

10 The {1 scale dependence in zq(^i) and 24(^4) induced by the Q2-Q4,6 mixing (multiplied by the factorized piece of 
/4 6 m Eq. 16.131 ) is cancelled by the fi scale dependence of in Eq. 16.141 . 
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Table 1: Numerical results for Reg 8 | „ „ for various input choices (see text). 



|{#)0i = 2 GeV)! 1 / 3 


L 5 = 1.0 x 10- 3 


L 5 = 1.5 x 10~ 3 


L 5 = 2.0 x 10~ 3 


0.260 GeV 
0.240 GeV 


0.35 - 1.23 
0.27- 1.24 


0.44- 1.21 
0.27- 1.23 


0.59 - 1.17 
0.41 - 1.22 



From these results, we conclude that a fair estimate of Reg 8 , at the 0(N C ) we have been working, lies 
in the range 

Regs = 1.33 ± 0.40 + 0.8 ± 0.4 = 2.1 ± 0.8 . (6.18) 

Q 2 ,Qi Q 2 -Eye 

In spite of the large errors involved, we find this result rather encouraging. When compared to the 
value |g 8 | oxp ~ 3.3 to be explained, it certainly points in the right direction towards a dynamical 
understanding of the AI = 1/2 enhancement. 



7 Phenomenology of e'/e 

In terms of the isospin amplitudes Aq and A 2 in Eqs. (|6.1[) and (|6.2(l . and to a very good approximation, 
one can write the CP violation observable e'/e as follows, 



£ 6 V2\e\ 



Inij42 ImAo 



Re^ 2 Re^o 



with 



r r , n j ReA 2 

$ = d 2 -oo + T — 0; and u = ^ a ■ 
4 Re^4o 



(7.1) 



(7.2) 



Using now the effective weak Hamiltonian in Eq. Ijl.3|l . one can obtain a formal expression of e'/e in 
terms of weak matrix elements of the four-quark operators Qi as follows 0] 



£ j=Im(V t * s V td ) 

e 2\e\\ReA \ 



p(°)(l-tt IB )--P (2) 



where 



P m =^Vi(p)((™)l\QiQj)\K°), for 1-0,2. 



(7.3) 
(7.4) 



and f^iB is a term induced by the effect of isospin breaking (m„ ^ rrid) 

1 (W 2 )ib 



IB 



(7.5) 



It turns out that, because of the AI = 1/2 enhancement factor 1/uj, in front of P^ on the r.h.s. 
of Eq. (|7.3() , and because of the values of the Wilson coefficients y$ [p) , the two P^ factors can be 
approximated, to a sufficiently good accuracy, as follows 



and 



P^ ^y s ^)((mr) 2 \Q 8 ^)\K ), 
where Qs denotes the electroweak penguin operator 

Qs = -12 ^2 e q^LqR){qii.d L ) . 

q=u,d,s 



(7.6) 
(7.7) 

(7.8) 
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with e q the electric charge of the quark q in e units. 

Recall that, since r in Eq. i|6.12[) is complex, the imaginary part of Cq and C4, and hence Img 8 |q 6j q 4 
is proportional to ye and 2/4. Our calculation of <?g| in sectional allows, therefore, for an 

evaluation of P' ' at the corresponding approximation i.e., lowest order in xPT and next -to-leading 
leading order in the 1/N C expansion, including terms of O(jf-), with the result 



P(°) ~ V2F {M% - ml) [y 6 {m c )f 6 (m c ;Mhad.) + Vi{rn c ) ferric ;/Mhad.)] 



(7.9) 



with the long-distance factors fe{m c ; /Zhad.) and fi{m c ; /-tiiad.) given in Eqs. 1)5. 2|l and (|5.3ll . 

We can also obtain an estimate of P^ from our previous work in ref. As discussed there, to 
lowest O(p ) in the chiral expansion, the four-quark operator Qs bosonizes as follows 



(Q 8 )|o(p») = -12 (O a (At)) tr (t/Af } C/tQ 



(7.10) 



<0|(ax,«ii)(dRdi)|0) 



where A^ = ^2^3 and = diag.[(2/3, — 1/3, — 1/3]. The vcv (<3 2 (a*)) also appears in the Wilson 
coefficient of the l/Q e term in the OPE of the Hlr(Q 2 ) correlation function 

l x e** x (0\T (u L ^d L {x)u Rl v d R {^) |0) = htftf - g^q 2 )U LR (Q 2 ) , (7.11) 

for which the MHA to large-iV c QCD gives a rather good approximation, as discussed e.g. in ref. |10| . 
This offers the possibility of obtaining an estimate of the vev 02(m) beyond the strict large- A^ c 
approximation, where O2 j('ipijj) 2 , and, therefore, without having to fix a value for the {^ip) 
condensate, which is poorly known at present. To lowest order in the chiral expansion (Q(p ) in this 
case as seen from Eq. I|7.10[) 'l we find [571155] 



P 



(2) 



, , Ml Ml 



25/8 
21/8 



NDR schen 
HV schem 



(7.12) 



We have now all the ingredients for a numerical evaluation of e'/e in Eq. J7J2J. For that purpose, we 
shall use the values of the physical parameters given in Appendix A. By far, the most sensitive param- 
eter in our determination of e'/e is the QCD quark-condensate. Low values ((^ip) 1 / 3 ~ —240 MeV) 
are favoured by various QCD sum rules determinations, like e.g. the determination in ref. |36| . using 
r-data. High values ((-tpip) 1 ^ 3 ~ —260 MeV) are favoured, however, by some of the lattice QCD 
simulations pJH] ■ For two recent determinations see, however, ref. [ID]. We, therefore, restrict the 
input value of (ipip) to a range 



|(^)(M = 2GeV)| 



1/3 



(250 ± 10) MeV. 



(7.13) 



Let us recall that the dependence of e'/e on (ijjip) appears in the term P(°\ trough the bosonization of 
the Qq operator i.e., the function f&(n; Hhad.) in Eq. 1)5.2(1 . and it is the sixth power of (-tpi/j) 1 ^ 3 which 
counts! 

Another important input parameter is the low-energy constant L5. As we have already discussed, 
the -L5 which appears in the factorized contribution induced by the Qe operator has to be taken as 
running, while L5 in the unfactorized contribution is constant, at the level of accuracy that we are 
working in the l/A c -expansion. As we have done in the calculation of Reggj^ Eye reported above, 

we leave i5|i arg c-jv c vary in the range 1 x 10 -3 < L 5 \ larso _ Na < 2 x 10~ 3 . 

Concerning the factor flis, we have nothing to add at present. We recall, however, that the 
estimate has changed from 

= 0.25 , See ref. gT] (7.14) 

to the more recent one xPT estimate [4*2"] 



rq B v = 0.16 ±0.03, 



(7.15) 
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which is the value that we shall be using here. We then have that 



- 7 = lm(V*V td ) 



G f uj 



2 e Re^o 



(0.055±0.008)xGoV~ 3 



pC) (i - n m ) - - p {2) 

0.84±0.03 S ~^T 



(7.16) 



We find large contributions for both p(°) , mostly induced by the Qe operator (the contribution from 
the Q4 operator amounts to less than 3% of the total), and for P^ 2 \ Typical numerical values are 
shown in Table 2 below for various input values of the quark condensate. The values in Table 2 have 
been obtained by letting i5|i argc - Nc and /^had. vary in the ranges 1 x 10 -3 < L^\ laiEQ _ Nc < 2 x 10~ 3 and 
0.8 GeV < /Xhad. < 1-3 GeV. The numbers correspond to the renormalization scale /i = m c = 1.3 GeV. 
The variations induced by the choice 0.8 GeV < [i < 2 GeV (in P^), or 0.7 GeV < M s < 1.1 GeV, 
change these numbers by 20%. 



Table 2: Numerical results for P(°'(l — fiie) in GeV 3 , (l/w)P^ 2 ' in GeV 3 and e'/e for various input 
choices of (ipip)- Here 1 x 10~ 3 < i 5 |ia Ig c-jv c < 2 x 10~ 3 and 0.8 GeV < ^ had . < 1.3 GeV. 



\{iM>)(ji = 2 GeV)! 1 / 3 


p(°i(i-n ]B ) x io 2 


Ip(2) x 10 2 


^ x 10 3 

e 


0.260 GeV 


8.1 - 9.2 


3.9 


2.0 - 3.4 


0.250 GeV 


6.3 - 7.2 


3.9 


1.2 - 2.1 


0.240 GeV 


5.0 - 5.6 


3.9 


0.5 - 1.1 



The results in Table 2 show that, for the range of values of the quark condensate (ipip) in Eq. I|7.13|l . 
the values we obtain for e'/e are perfectly compatible with the latest world average 03] from the 
NA31, NA48 and KTeV experiments 

Re(e'/e)| Exp = (1.66 ± 0.16) x 10~ 3 . (7.17) 

A vanishing or negative value of e'/e, as reported by the lattice QCD groups |441 145| . appears to us 
very unlikely. 

For the purpose of comparison with other theoretical predictions, it is convenient to represent 
Eq. H7.16J) as a straight line in a plane with plotted in the horizontal axis, and P(°)(l - n IB ) 

in the vertical axis This is the plot shown in Fig. 8, where the width of the solid line reflects 

the experimental errors in the overall factor in the r.h.s. of Eq. I|7.16(l and in Eq. I|7.17|l . The cross in 
the same plot, represents the prediction which follows from our present estimates of P^ and P^ 11 , 
for the restricted input of the quark condensate in Eq. (|7.13(l . but including the errors of scanning 
the other input parameters in Eqs. (A. 10) to (A. 11) of the Appendix, and varying the matching scale 
between 800 MeV and 2 GeV. Furthermore we have also allowed for violations of factorization in the 
residue of the OPE in Eq. (|4.20|1 by a factor of 2, as we have done before for Eq. H4.21|l . It turns 
out, however, that the influence of a possible error in this residue is now much milder because the UV 
behaviour of the integral of Q 2 Wdgrr{Q 2 ) is very much dominated by the Goldstone double pole 
term. 



8 Discussion and conclusion 

We hope to have shown how one can treat analytically, and on a sound theoretical basis, the calculation 
of electroweak matrix elements beyond factorization. In particular, one important goal of this work 

11 The error in p( 2 ) does not include other analytic determinations 1491 15UI lo 1 1 1521 1531 and lattice QCD determina- 
tions I54ll55ll44ll45) . 
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was to study the magnitude of the unfactorized contributions in two cases of special interest, i.e. e'/e 
and the AI =1/2 rule. In section 4, we explained how, bringing new hadronic scales, the unfactorized 
contributions may turn out to be sizeable as compared to the factorized ones, depending on the scales 
involved. Already at the level of the leading 0(n,f/N c ) corrections, we have shown by an explicit 
calculation that this is indeed the case for the Qe and Q4 penguin operators. 

A 



12 
10 
8 

6 
4 

? 



1 2 3 4 5 



Fig. 8 The experimental value of e'/e in Eq. [7.17}) fixes the straight line in the the figure, with a 
certain width due to experimental errors, see Eq. [7.16}) . Our theoretical prediction, including errors, 
is represented by the cross with error bars. 

For the operator Qq, the unfactorized contribution turns out to be even much larger than its 
factorized one, due to the large chiral coefficient ~ L5 — in Eq. (|4.18() which through Eq. (|4.29() 
brings a large coefficient to the hadronic scales in Eq. (|5.2|1 . As a result, the leading piece is not 
the factorized contribution but the unfactorized one for which we just calculated the lowest order 
contribution. From this leading lowest order contribution, the higher order corrections should now be 
calculated in subsequent works. Similarly for Q4 (and for the directly related eye contribution of Q2), 
we find a rather large unfactorized contribution through Eqs. (|4.30l) and l|5.3|l . 

Note that the claim that unfactorized contributions are large was already made in refs. |17l 1191 
EH EH for Qi, 2 and in refs. 0123122] for Qe- In refs. [1112131, the calculat ion of the unfactorized 
contribution was done for Qq by just taking the effects of the pseudoscalar mesons regularizing the 
pseudoscalar loops with an euclidean cut-off A c . From our present analysis, which incorporates the 
effect of hadronic resonances explicitly, we observe that this can be justified at the qualitative level, 
because, as was anticipated in those references, a quadratically divergent ( ~ A^) piece can be inter- 
preted as the reflect of terms proportional to the square of the masses of the hadronic resonances. In 
our more refined approach, the cut-off K 2 C is replaced by combinations of M v , M|, Mp, and /i^ad ■ 
This is also what happens for the quadratic coefficient of the Qq matrix element (i.e. the r.h.s of 
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Eq. 1)4. 29[) proportional to Lg — ^L^). 12 However, there is no way, with the cut-off procedure advo- 
cated in refs. |171ll8lH51l20| . to fix unambiguously the choice of A c , and there is no control at all on 
the short-distance matching scale dependence. The unfactorized contribution has also been estimated 
in refs. 0Jj within the constituent chiral quark model of ref. \'2'A\ . again without short-distance match- 
ing. In refs. |21l I22| . the previous problem is somewhat ameliorated, although within the framework 
of the ENJL model. It turns out, however, that in this model it is also found that the unfactorized 
contributions are large. 

The advantage of the large--/V c QCD approach presented here is that it allows to use the same 
dimensional regularization both at short and long distances. We have been able to do a fully analytic 
matching of scale dependences, which allows for a clear separation of the various scales and hadronic 
quantum numbers involved in the problem. 

To improve our calculation requires the incorporation of several effects which are beyond the 
scope of the present analysis. The 0(l/N c ) not enhanced by a nf factor should be calculated to 
have a full NLO calculation in the l/N c expansion. The long-distance scheme dependence should 
also be calculated to cancel the corresponding short-distance one (which can be done exactly in our 
framework as was shown for Bk in ref. m]). The hadronic ansatz could be refined, systematically, 
by considering higher dimensional short- and long- distance constraints. Moreover, the effect of the 
final state interactions (FSI) could eventually be incorporated, following the work in ref. |48|. Notice, 
however, that the effect of the unfactorized contributions which we find is much larger than the size 
of the FSI effects estimated in this reference. 13 

Finally, we observe that a precise calculation of e' /e, within the framework we have discussed 
here, is correlated with the determination of some of the low energy constants entering our input; 
in particular the value of the quark condensate (ipip) and to a lesser extent the L5 coupling. Other 
important input values are f2/s (where it would be nice to include as well the unfactorized contribution 
which, so far, has been always neglected), Aqcd and Im Vt s V t * d . 

As a result of all the considerations discussed above, we believe that a systematical analytic 
calculation of e'/e in the Standard Model is now becoming conceivable although, in our opinion, 
matching the level of the experimental precision is likely to become a very difficult task. 
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APPENDIX 

A Compilation of Numerical Inputs 

The numerical input we have used in the text are as follows 14 : 

Im (V t * s V td ) fUJ = 0.055 ± 0.008 , Ref.[3S]. (A.l) 



2e 



Re^o 



12 Note however that as was already pointed out above the important L3 contribution was missing in this quadratic 
coefficient in ref. 1201 . As a result the enhancement of the Qa contribution to e'/e found in ref. 1201 was smaller than 
the one we obtain here. 

13 In refs. [481 , the unfactorized contribution accounts for an error of ±0.5 • 10 — 3 in their quoted value of e'/e. We find 
here that this is an underestimate. 

14 Notice that the Aq amplitude in ref. |3] is ours times a factor W3/2. 
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and 



M w = 80.4 GeV, sin 2 6 W = 0.23 , a = 1/137.036, (A.2) 

m c = 1.3GeV, m h = 4.4GeV, m t = 170GeV, (A.3) 

Mk — 0.498 GeV , = 0.135 GeV , u = 1/22.2, (A.4) 

i^o = 0.087 GeV, My = 0.770 GeV , Ms = 0.9 ± 0.2 GeV , (A.5) 

M A = 1.2 GeV, M v > = 1.4 GeV , M F / = 1.3GeV, (A.6) 

L 3 = -3.0 • 10~ 3 , L 5 (M P ) = 1.4- HP 3 , L 9 (M P ) = 6.9 • 10 -3 . (A.7) 



(V>V>)(At = 2 GeV) = -(250 ± 10 MeV) 3 , (A.8) 
i 5 |la rB e-« c = (1.5 ± 0.5) x 10~ 3 (A.9) 



MS 

= 0.16 ±0.03. (A.ll) 



Ai^L = (325 ±40) MeV, (A.10) 



The Wilson coefficients we have used have been calculated using the same program as the one which 
has been used in ref. 0). For all calculations we used LO Wilson coefficient except for the evaluation 
of P( 2 ) in Eq. $77^1 for which we took the NLO-NDR Wilson coefficients. 
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